Local structure evolution in polycrystalline Zni.^jMg^O (0 < x < 0.15) studied by 
Raman and by synchrotron x-ray pair distribution analysis 
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The local structures of Zni_j,MgiO alloys have been studied by Raman spectroscopy and by 
synchrotron x-ray pair distribution function (PDF) analysis. Within the solid solution range (0 < 
X < 0.15) of Zni-^MgiO, the wurtzite framework is maintained with Mg homogeneously distributed 
throughout the wurtzite lattice. The iSj'^'^ Raman line of Zni_^Mga;0 displays systematic changes 
in response to the evolution of the crystal lattice upon the Mg-substitution. The red-shift and 
broadening of the iJj'^'^ mode are explained by the expansion of hexagonal afc- dimensions, and 
compositional disorder of Zn/Mg, respectively. Synchrotron x-ray PDF analyses of Zni-ajMgxO 
reveal that the Mg atoms have a slightly reduced wurtzite parameter u and more regular tetrahedral 
bond distances than the Zn atoms. For both Zn and Mg, the internal tetrahedral geometries are 
independent of the alloy composition. 

PACS numbers: 61.10.Nz, 71.55.Gs, 77.22. Ej 



I. INTRODUCTION 

Polar semiconductors such as wurtzite ZnO and GaN 
have great potential for use in the polarization-doped 
field effect transistors (PolFETs) which exploit a polar- 
ization gradient at the channel layer to attain higher mo- 
bility and higher concentrations of carriers. The per- 
formance of PolFETs depends primarily on the interface 
quality and the magnitude of the polarization gradient at 
the heterojunction. Those goals can be approached by in- 
terfacing a polar semiconductor with its alloy derivative, 
as exemplified by the fabrication of a PolFET compris- 
ing the interface of GaN/Gai_^AUN (a; = - 0.3). ^ In 
designing ZnO-based PolFETs, Zni_a;Mga;0 alloys can 
be considered as the sub-layer component, since Mg- 
substitution effectively alters the polarization of ZnO 
while keeping the lattice dimensions nearly unchangedi^ 

In a previous report,^ we have outlined the prospects 
for ZnO/Zni„a;Mg2;0 heterojunctions for PolFET appli- 
cations, based on structural analysis of Zni_a;Mga;0 al- 
loys using synchrotron x-ray diffraction. We showed that 
the ionic polarization can be tuned by «14% from ZnO to 
Zno.85Mgo.15O, despite only small changes in the cell vol- 
ume (wO.3%). Also, it was found that Mg can substitute 
up to 15% of total Zn in the wurtzite lattice, without any 
evidence for segregation. While the average crystal struc- 
tures of Zni_a;Mga;0 phases were accurately determined 
in the above study, there still remain the local structural 
details of the solid solutions to be understood. Therefore 
in a continuing effort, we have performed Raman studies 
and synchrotron x-ray pair distribution function (PDF) 
analysis of polycrystalline Zni_2,Mg2,0 samples. Raman 
spectroscopy is known to be useful for examining com- 
positional disorder and/or the presence of strain within 
semiconductor alloys i^i^i^ii^ In previous Raman studies 
on wurtzite type crystals, the peak position and shape of 
the E^^^^ phonon mode have been used to investigate the 
effects of sample grain size, heterogeneous components. 



and defects. ^'^ PDF analysis provides a powerful probe 
of local non-periodic atomic displacements in the short 
range, and is a suitable complement to the /c-space re- 
finement technique. ^'^^ This ability is very relevant to 
the possible distinction of tetrahedral out-of center dis- 
placements of Zn and Mg, and provides for a better com- 
parison platform with density functional calculations of 
the crystal structure. 

In this study, we provide a complete description of 
polycrystalline Zni_a;Mga;0 alloys prepared from crys- 
taUine Zni_2;Mga;(C204)-2Il20 precursors. We present 
results of thermogravimetry, ultraviolet/ visible diffuse- 
reflectance spectroscopy, infrared, and Raman spec- 
troscopy in addition to examining the distinct tetrahe- 
dral geometries of Zn and Mg using synchrotron x-ray 
PDF analyses based on supercell structure models. 



II. EXPERIMENTAL 

Powder samples of Zni^JAgr^O [x = 0, 0.05, 0.10, 
and 0.15) were prepared using an oxalate precursor route 
as described previouslyi^ Aqueous solutions of zinc ac- 
etate, magnesium acetate, and oxalic acid were sepa- 
rately prepared and mixed to precipitate zinc magnesium 
oxalates, which are crystalline, single-phase compounds 
with Zn^+ and Mg^+ homogeneously mixed at the atomic 
level. The precipitates were thoroughly washed with 
deionized water and dried at 60° C for 4h to produce 
white powders of Zni_a;Mga:(C204)-2Il20, as verified by 
powder x-ray diffraction. Subsequently the oxalate di- 
hydrates were transformed to 'Z^ni—x MgxO by heating 
in air at 550°C for 24 h. The thermal decomposition 
of Zni_2,Mg2;(C204) -21120 was monitored by thermo- 
gravimetry using a Cahn ThermMax 400 thermogravi- 
metric analyzer (Thermo Scientific). For each composi- 
tion, «60mg of powder was heated in air up to 1000°C 
at 5°Cmin"^. Fourier-transform infrared (FT-IR) spec- 
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tra of Zni_2;Mga;0 powders were recorded in KBr using 
a Nicolet Magna 850 FT-IR spectrophotometer in the 
transmission mode. Diffuse-reflectance absorption spec- 
tra were measured for Zni_a;Mga;0 in the wavelength 
range of 220— 800 nm using a Shimadzu UV-3600 spec- 
trophotometer equipped with an ISR-3100 integrating 
sphere. The powder samples were mounted to have 
flat surfaces and «1.5mm thickness. The optical band 
gap was determined by extrapolating the linear part of 
absorption edge to zero-absorption level. Raman mea- 
surements were conducted at room temperature using 
an optical microprobe fltted with a single monochroma- 
tor (Jobin-Yvon, T64000). Powder samples were com- 
pacted on frosted glass plates and spectra were recorded 
in backscattering geometry using a 488.08 nm Ar"*" laser 
with a beam power of 50 mW and a spot size of w2 /im. 
For each sample, five acquisitions of 30 s exposure were 
performed and the averaged spectra are reported. The 
E'a'^'^ phonon mode was chosen for detailed peak profile 
analyses. The spectral background was removed follow- 
ing Shirley^ii and Breit-Wigner type peak fitting^^ was 
employed to determine peak position and the width. 




FIG. 1: Thermogravimetry of Zni_^Mg:,0 {x = 0, 0.05, 0.10, 
and 0.15) in air at a heating rate 5°Cmin~^. 



III. RESULTS 



Synchrotron x-ray scattering experiments were carried 
out at beam line 1 1-ID-B of the Advanced Photon Source 
(Argonne National Laboratory) using an x-ray energy of 
90.8 keV (A « 0.1365 A) at room temperature. Use of 
high energy radiation enables data acquisition at high Q 
(=47rsin0/A) wave vector, which in turn improves the re- 
liability of the Fourier transformation for obtaining the 
PDF G{r). In the present work, scattering data with 
Qmax of 28 A^^ were utilized for extracting PDFs. Sam- 
ple powders were loaded in Kapton tubes and data were 
measured in the transmission mode using an amorphous 
silicon image plate system (General Electric Healthcare). 
For each sample, 33 images were taken with an exposure 
time of 16 s per image. The program FiT2Di^ was used to 
convert images into the corresponding one-dimensional 
x-ray diffraction (XRD) pattern. The average crystal 
structures of these samples have been determined by Ri- 
etveld method and reported previously)^ For obtaining 
the PDF from XRD data, the program PDFgetX2 was 
usedM- First the measured scattering intensities were cor- 
rected for sample container background, Conipton scat- 
tering, and Laue diffuse scattering. Next the coher- 
ently scattered intensities I{Q) were normalized in ab- 
solute electron units to give total scattering structure 
functions S{Q). Finally the reduced structure functions 
F{Q) — Q[S{Q) — 1] were Fourier-transformed to produce 
the atomic PDFs G(r). The refinements of Zni-x^gxO 
structures were performed against the above obtained ex- 
perimental PDFs using the software POFfiTf^^ 



A. Sample characterization 

The thermogravimetric profiles of 

Zni_^Mg^(C204)-2H20 {x = 0, 0.05, 0.10, and 
0.15) are shown in Fig.[T] The stepwise processes of 
dehydration and oxalate-to-oxide conversion are clearly 
observed at around 130 and 390°C, respectively^^ with 
the weight changes in good agreements with calculated 
estimates. With increase of Mg content, the onset 
temperatures of both weight loss steps shifted to higher 
temperature. This could be interpreted as indicative 
of the extra activation barrier for stabilizing Mg in 
the wurtzite lattice. For all cases, oxide formation 
was nearly complete at 500°C without any noticeable 
weight change at higher temperature. It is therefore 
assumed that isothermal heating at 550°C for 24 h 
completely transforms Zni_a;Mga;(C204)-2Il20 to the 
oxides Zni„a;Mga;0. 

In Fig.d the FT-IR spectra for Zni_a;Mg^O are 
shown, which are consistent with the registered reference 
data of wurtzite ZnOf^i A previous report suggested that 
organic precursors for ZnO synthesis may leave carbon- 
ate species strongly bound within the lattice,^* but the 
samples in this study did not exhibit any spectral fea- 
ture at wl300 and wl500cm~^ demonstrating the com- 
plete combustion of oxalate,— The wurtzite lattice vibra- 
tions were observed as broad IR bands at 400— 600cm~^. 
Upon Mg-substitution, these stretching modes shifted to 
higher wavenumber as a result of the smaller reduced 
mass of Mg— O compared with that of Zn— O. 

The optical band gaps (Eg) of Zni_a;Mga;0 powders 
were determined from the diffuse-reflectance absorption 
spectra as shown in Fig. [31 The absorbance (A) of the 
samples were obtained from the measured reflectance (i?) 
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FIG. 2: FT-IR spectra of Zni_^Mg:,0 {x = 0, 0.05, 0.10, and 
0.15). 
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FIG. 3: Diffuse-reflectance absorption spectra for 
Zni_j;Mg:rO {x = 0, 0.05, 0.10, and 0.15). For each spectrum, 
the band gap energy is indicated by the extrapolated arrow. 

according to the Kubclka-Munk relation.— 



In the present measurement setup, the sample analytes 
are thick enough («1.5 mm) to disallow transmission, and 
the Kubelka-Munk theory is appropriate. As clearly seen 
in Fig. [31 the absorption edges are shifted to the higher 
energy side with the increase of Mg content x. Also 
noticed are slight broadenings of the edge slopes upon 
the Mg-substitution. The band gap energy gradually in- 
creased from 3.24 {x = 0) to 3.26 {x ^ 0.05), 3.30 {x 
= 0.10), and 3.35 eV {x = 0.15). In previous studies 
on Zni_a.Mga;0 thin films grown by pulsed laser deposi- 
tion, Ohtomo etal. have observed a similar composition- 
dependence of the band gaps; monotonic increases from 
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FIG. 4: Raman spectra of the Zni_^Mg^O {x = 0, 0.05, 
0.10, and 0.15). A plasma line from the 488.08 nm Ar^ laser 
appears at Ri348cm^^ (indicated by a vertical line). 

3.30 (x = 0) to 3.63 eV (x = 0.14)^ Their absolute band 
gap energies are not exactly reproduced in our polycrys- 
talline samples, but this may be due to the differences in 
experimental details. The optical band gaps of solids are 
frequently determined to be different depending on sam- 
ple type, optical characterization technique, and band 
gap determination methodi ^^'^'^ In the above study on 
Zni_a;Mgj;0 films the band gaps were evaluated with 
the assumption oc {hv — Eg), in which the absorp- 
tion coefficient a was deduced from the transmittance 
measurements. 



B. Raman spectroscopy 

Wurtzite ZnO belongs to the space group 
(Hermann-Mauguin symbol PQ^mc) and has six lattice 
phonon modes, Ai + 2Bi + Ei + 2E2, of which the Bi 
branches are Raman inactivc i^^i^^ The Ai and Ei modes 
are polarized along the z-direction and xy-plane, respec- 
tively, whereas the two E2 modes (£^2°^ and £^2'^'^) 
non-polar Each of the polar modes is split to lon- 
gitudinal (LO) and transverse optical (TO) components 
due to the macroscopic electric field associated with the 
LO phonons. 

Raman spectra of polycrystalline Zni_2;Mga;0 powders 
were examined over the frequency range 90— 830cm~^ 
and part of those spectra are presented in Fig.|4l For all 
cases ofx = — 0.15, the measured Raman spectra agree 
well with the wurtzite ZnO vibration modes, without any 
new bands arising from Mg-substitution. Raman features 
characteristic of cubic MgO^ were not observed. In pre- 
vious Raman studies on ZnO films or nanophases doped 
with N, AI, Ga, Sb, Fe, Mn, or Mg, additional bands have 
been observed and attributed to the induced lattice de- 
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FIG. 5: (a) phonon frequency cuo and (b) width F of Raman 
-Ej'^ mode as functions of Mg content in Zni-^^MgajO. Gray 
lines are best linear fits. 



feet or the dopants' local vibration . A previous Ra- 
man study of polycrystallinc Zni_a;Mga;0 has shown that 
the Mg-substitution flattens the overall Raman signal. 
In Fig.H the Raman peaks from Ai(TO), Ei{TO), E^°^ 
(multiphonon), and E^^'^''^ modes are found. The Ai and 
El modes reflect the strength of the polar lattice bonds, 
which are of interest in relation with the c-axial displace- 
ments of cations. However the peaks corresponding to 
those two modes were not well resolved from the back- 
ground and could not be used for a detailed study. Usu- 
ally Raman signals from polar modes are weaker in in- 
tensity due to phonon-plasmon interactions!^ The E^^^^ 
phonon, which is the most prominent in the ZnO Ra- 
man spectra, was used for the peak profile analysis. The 
Raman line shape was fitted using the Breit-Wigner func- 
tion. 



[1 + 2P{lu - ujo)/TY^ 
1 + [2{lj - c^o)/r]2 



(2) 



where I{uj) is the Raman intensity at a given frequency, 
ujQ is the Raman shift, F is the broadening expressed as 
full- width at half-maximum (FWHM) , and f3 is an asym- 
metry parameter (/3 = for the symmetrical peak):l- For 
each Raman spectrum of Zni-^^Mg^:©, the £^2'*'' phonon 
frequency loq was determined along with F, as plotted in 
Fig.El 

If we consider the mass change in the Mg-substitution 
for Zn, the phonon frequency is expected to increase for 



the samples with higher x. But the observed trend is 
contrary to this: The E'j'®'' phonon frequency decreases 
upon the Mg-substitution. To explain this, it should be 
recalled that the E2^^^ mode of ZnO corresponds mostly 
(«85%) to the vibration of oxygen atoms^^ and is in- 
sensitive to the mass substitution on the cation site. 
On the other hand, since the iJj'^'^ mode is associated 
with atomic motions on the xy-plane, its phonon en- 
ergy depends on the in-plane lattice dimensions. In the 
Zni_a;Mgj:0 solid solution, the hexagonal parameter a 
monotonically increases with x, hence softening the E2 
phonon mode. The decrease of £'2 '^'^ mode frequency has 
been similarly observed from ZnO thin films under the 
tensile strain along the xy-planc.'^'^ 

In addition to the frequency shift, it is also noted that 
the £2'^'^ peaks are markedly broadened upon the Mg- 
substitution. The Raman mode frequency and the peak 
width can be influenced by the random substitution of Zn 
with Mg through the compositional disorder effect i^i^i^^ 
As inferred from the evolution of average crystal struc- 
tures in Zni_a;Mga;0, the Mg atoms do not have an iden- 
tical bonding geometry to that of Zn. In other words the 
Mg-substitution alters the translational symmetry of the 
wurtzite lattice to modify the phonon oscillation fleld. In 
previous Raman studies on the alloy systems, the spatial 
correlation has been regarded an adequate mechanism 
responsible for the red-shift of the mode frequency as 
well as the peak broadeningjSiSi^^ The model is based on 
the wave vector uncertainty, Aq — 271 /L, where q is the 
phonon wave vector range, L is the length scale of phonon 
confinement, and q has the phonon dispersion relation 
Lo{q) with the Raman frequency. It predicts that, as the 
L becomes smaller, Aq becomes larger so that a wider 
range of frequencies is allowed for Raman scattering. 

In Zni_2;Mga;0 compounds as well, the presence of 
the hetero component Mg seems to perturb the wurtzite 
phonon modes. As shown in Fig.[5l both the frequency 
shift and the peak broadening become greater for the 
samples with higher Mg-concentration. However, the 
symmetry of i?2 '^'^ peak shape did not show significant 
dependences on the Mg-concentration, which is some- 
what unexpected. This is in contrast to the spatial cor- 
relation model and also with the experimental examples 
of disordered semiconductor crystals, where the Raman 
line shapes become more asymmetric as the peak widths 



increase 
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C. Pair distribution functions 

Both Rietveld refinement and PDF analysis utilize one- 
dimensional diffraction patterns to extract crystal struc- 
ture information. The former focuses on the long range 
average structure based on Bragg diffraction, while the 
latter focuses on the local structure information over 
short length scales from the diffuse scattering in addition 
to the Bragg diffraction^ii^i^i^ii^ Therefore the above 
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FIG. 6: XRD patterns of Zno.g5Mgo.05O, used for (a) PDF 
analysis and (b) Rietveld refinement. Sample-to-detector dis- 
tances were 150 mm for (a) and 660 mm for (b). 

two methods rely on distinct data qualities. The dif- 
fuse scattering becomes more prominent at high Q, for 
which the PDF data require very large Ewald spheres, 
but not necessarily very high resolution, and extremeley 
large signal to noise ratios. For Rietveld analyses, highly 
resolved Bragg peaks are preferred, at the expense of col- 
lecting data over smaller Ewald spheres. In order to meet 
both needs, the diffraction experiments were carried out 
at two different sample-to-detector distances, 150 mm for 
the PDF data and 660 mm for the Rietveld data. Fig. [5] 
compares the two XRD patterns of Zno.95Mgo.05O, one 
used for Rietveld refinement and the other for PDF study. 
The pattern for Rietveld refinement comprised sharp and 
well-resolved Bragg peaks but the diffraction information 
rapidly vanishes at Q beyond 18 A"-'^. However, in the 
pattern for PDF analysis the scattering intensities are 
observable at Q as high as 30 A"-'^. 

The PDF data G(r) for the four samples of 
Zni_::cMga;0 {x = 0, 0.05, 0.10, and 0.15) were obtained 
following a standard data processing sequenceiJ--'-^ In 
Fig. [3 the PDFs from different Zni_^Mga;0 compositions 
are superimposed, and it is observed that the data are 
hardly distinguished. Since the average crystal structures 
(lattice symmetry, lattice parameters, atomic positions, 
etc.) of Zni_2;Mg2,0 phases are similar, the composition- 
dependent changes in PDFs are not immediately notice- 
able. However a close view at the PDF peak at r w 2 A 
exposes a subtle evolution in the coordination shell of 
the nearest cation— anion bond pairs. As the Mg content 
increases, these PDF peaks are slightly shifted to longer 
r, with gradual increases of peak heights. This finding 
is consistent with the previous Rietveld refinement result 
which indicated that the (Zn,Mg)— O bond distances be- 
come more regular upon Mg-substitution.^ The first shell 
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FIG. 7: Experimental PDFs for Zni_a;Mga;0 phases with x 
— 0, 0.05, 0.10, and 0.15. Inset shows the zoomed-in view of 
the peak corresponding to the first coordination (Zn,Mg)— O 
shell. 



interatomic distances are (3x)1.973lA and 1.9941 A in 
ZnO, and (3x)1.9756A and 1.9873A in Zno.85Mgo.15O. 
Although small in magnitudes, the above changes corre- 
spond to «45% reduction in the standard deviation of the 
four bond distances, which result in the discernible peak 
sharpening. The shifts of peak positions can be ascribed 
to the lengthening of the three shorter bonds rather than 
a change in the longer apical bond. 

The real space PDF data are therefore useful for di- 
rectly depicting the bonding geometry in a specific r- 
range. For more detailed and quantitative information, 
the PDF data were analyzed by full profile refinements 
using several different models of simple cell as well as 
supercell. The simple cell models were based on the 
hexagonal wurtzite structure (space group PQ^mc, Z = 
2, a « 3.25 A, c « 5.20 A) with O at (i §0) and Zn/Mg 
at (i|u). In the analyses using the simple cell, the 
fractional occupancies of Zn and Mg were set accord- 
ing to the composition of each phase Zni_j,Mgj,0 {x — 
0, 0.05, 0.10, and 0.15). Supercell model structures were 
constructed by expanding the wurtzite simple cell to a 
3x3x1 size (Z = 18, a 9.75 A, c « 5.20 A) . By partially 
replacing the Zn atoms with Mg, the parent ZnigOis su- 
percell can be modified to the compositions ZniyMgOig 
(Zno.944Mgo.056O), Zni6Mg20i8 (Zno.889Mgo.111O), and 
Zni5Mg30i8 (Zno.833Mgo.i67O) to model the x = 0.05, 
0.10, and 0.15 phases, respectively. The Mg atoms in 
Zni6Mg20i8 and Zni5Mg30i8 models were separated to 
be as distant from one another as allowed by the cell. 
PDF refinementsi^i^ of G(r) were carried out using the 
variables of scale factor, dynamic correlation factor, res- 
olution factor, lattice constants, independent u of Zn and 
Mg, and isotropic temperature factors of Zn, Mg, and O. 

First the PDF refinements were performed for 
Zni_3;Mg^0 {x = 0, 0.05, 0.10, and 0.15) using the four- 
atom simple cell models, in which the cation composi- 
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FIG. 8: PDF refinement profiles for (a) ZnO and (b) 
Zno.85Mgo.15O, using four-atom simple cell models. Exper- 
imental (open circles) and calculated (solid lines) data are 
plotted in the bottom panels, along with the difference pat- 
terns (Gcxpt — Gcaic) in the top panels. 
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FIG. 9: Position parameters for Mg and Zn in Zni_iMga;0 
phases, as analyzed by PDF refinements (1.5— 10 A) of simple 
cell and 3x3x1 superceil model structures. Rietveld refine- 
ment results^ are also compared. For the PDF results, error 
bars indicate one estimated-standard-deviation (esd) while for 
the Reietveld results, the symbol height corresponds to one 
esd. 



TABLE I: Summary of PDF refinements for Zui-^Mg^^O {x 
= 0, 0.05, 0.10, and 0.15) using simple cell models, over the 
r-range of 1.5— 10 A. 



X 





0.05 


0.10 


0.15 


iiw (%) 


13.4 


17.0 


17.2 


17.8 


a (A) 


3.2500(2) 


3.25028(5) 


3.2504(1) 


3.2506(1) 


c(A) 


5.2101(4) 


5.2109(1) 


5.2115(3) 


5.2119(3) 


u{Mg) 




0.380(4) 


0.379(4) 


0.380(3) 


u{Zn) 


0.3819(4) 


0.3817(1) 


0.3815(3) 


0.3813(3) 



tions were accounted for with fixed occupancy param- 
eters. TableU] summarizes the results of the 1.5— 10 A 
range fitting, and Fig. [5] shows the refinement profiles 
for ZnO and Zno.85Mgo.15O. Refinements over wider r- 
ranges tended to yield higher factors but the fit re- 
sults were similar to those given in TableHl It can be 
mentioned that the obtained i?w's near 20% are common 
for PDF refinements, even for well-crystallized materials. 
The i?w's from the PDF and Rietveld refinements cannot 
be directly compared, but they serve the same purpose 
of finding out the best structure solution from a number 
of competing models with similar numbers of refinable 
parametersi^2ii£ The inherently higher R„ in the PDF 
analysis stems from the greater sensitivity to the local 
atomic ordering, imperfect data correction, and system- 
atic errorsi^ The i?w's for the Zni_j;Mga;0 {x > 0) sam- 
ples are significantly higher than that for ZnO (cc = 0). 
It implies the presence of irregular Mg/Zn distribution 
in the alloy phases, which may not be well portrayed by 
simple structure models. 

For all the four compositions, the lattice constants 
from the PDF agree well with the Rietveld refinement 
results, although the precision is slightly lower in PDF. 
Moreover the Zn atomic position and the Zn— O bond dis- 
tances in ZnO are well reproduced in the PDF analysis, 
when compared with Rietveld refinement results. For the 
Zni_a;Mg2:0 phases, the z-coordinates of Zn and Mg were 
independently refined to investigate the distinct local ge- 
ometries of the two cation types. Interestingly the refined 
u{Mg) is smaller than M(Zn) for all the three solid solu- 
tions with X = 0.05, 0.10, and 0.15 (Fig.©. This trend 
is consistently repeated in the superceil model analyses 
detailed below. The u(Mg) parameters have rather large 
estimated deviations together with weak variations de- 
pending on X, but the observation u(Mg) < u(Zn) was 
quite robust regardless of the refinement r-range or the 
choice of constraints. 

We recognize that in the above simple cell model for 
Zni_a;Mgj;0, only the cation sub-lattices can be distorted 
while the anion framework remains undistorted. The 
wurtzite structure is its own antitype where the cation 
and anion sublattices are interchangeable. Depending 
on which of the two is distorted the resulting PDFs will 
be slightly different from each other, although the corre- 
sponding average crystal structures are practically indis- 
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(a) cation-lattice deformation 



(b) anion-lattice deformation 



FIG. 10; (Color online) Local structural distortions in Zni_a;Mg2,0 supercells of (a) cation-lattice deformation, M(Mg) 7^ u(Zn), 
w(OMg) = ii(Ozn) and (b) anion-lattice deformation, u{Mg) = u(Zn), M(OMg) / u(Ozn). ax and eq denote the cation— anion 
distances with different bond orientations, || c and _L c, respectively. As resulted from local distortions, ax' 7^ ax and eq' 7^ eq. 



tinguishable. In case the oxygen-lattice is distorted, the 
partial PDF from 0—0 pairs is modified from that of 
the average crystal structure, and likewise, distortion of 
the cation-lattice affects the PDF for (Zn,Mg)-(Zn,Mg) 
pairs. The differentiations of cation— anion bond lengths 
in Zni_a;Mga;0, resulting from those two deformation 
cases, are illustrated in Fig.[Tni In order to verify the 
lattice deformation type, we have examined the PDF re- 
finements of Zno.85Mgo.15O using two different supercell 
models of ZnisMgaOig. In the anion-lattice deforma- 
tion model, the ^-coordinates of Zn and Mg atoms were 
fixed to or i, but the z-coordinates of O atoms were 
refined with constraints; the three O atoms (OMg) that 
are apical to Mg were bound to have symmetry-related 



z-coordinates {zj 



i), and so were the other 



fifteen O atoms (Ozn)- In the cation-lattice deformation 
model, the z-coordinates of O atoms were fixed, while 
the Zn and Mg positions were allowed to vary in inde- 
pendent groups. For each model, PDF refinements were 
carried out for the r-ranges from 1.5 A to various r^ax 
of 6.4— 21. 4 A. Attempts using the smaller rmax (< 6 A) 
suffered instability of the lattice constants, and were dis- 
carded. Since the resolution factor was very sensitive to 



the 



t, it was fixed to 0.1 in all the refinements. 



The refinement results using the two supercell mod- 
els are presented in Fig.[TlJ For both models, the i?w 
factors were obtained in the range of 17'^19.5% with no 
apparent Tmax-dependences. The two models are based 
on the deformation of different sublattices, and accord- 
ingly the z-coordinates of either cation or anion are pre- 
sented for each. However, it can be well assumed that 
M(Mg) and u(Zn) from the cation-lattice deformation 
model are comparable respectively to u(OMg) and u(Ozn) 
from the anion-lattice deformation. As shown in Fig.fTTl 
the two supercell models similarly indicate that u(Mg) 
< u(Zn) and u(OMg) < u{Ozn), namely that the axial 
Mg— O bonds are shorter than the axial Zn— O bonds. 
A quick examination of Fig. [11] reveals that the cation- 
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FIG. 11: Atomic positions in Zno.85Mgo.15O as functions of 
the length of the vector rmax- PDF refinements were carried 
out using the supercell model structures based on (a) anion- 
lattice deformation and (b) cation-lattice deformation. Gray 
lines are best linear fits. 



lattice deformation model provides more reproducible re- 
finements of atomic positions. The averages over the 
''max-dependent refinement results are M(Mg) = 0.3795(6) 
vs. w(OMg) = 0.3765(10) and M(Zn) = 0.3813(1) vs. 
w(Ozn) = 0.3813(5). The comparison of lattice constants 



8 



o Zn-03,(x1) 
□ Zn-0^^(x3) 



M9-0,Jx1) 



1.99 - 



- 1.98 



o 



n □ 



■j I I I I I I I I I I I I I I 

5 10 15 

FIG. 12: Zn-O and Mg-0 bond distances in Zno.85Mgo.15O, 
obtained from the PDF refinements using the cation-lattice 
deformation supercell model. Gray lines are guides to the 
eye. 



(not shown) further showed that the anion-lattice defor- 
mation model yields more scattered results. It is there- 
fore judged that the cation-lattice deformation is appro- 
priate for analyzing the Zui-x^gxO structure. 

Here it can be instructive to compare the local dis- 
tortions in both deformation models in detail. In the 
cation-lattice deformation model (Fig.fTOk). the bond- 
ing geometries of cations are differentiated simply by the 
atom type, Zn or Mg. On the other hand the anion- 
lattice deformation (Fig.fTOb) creates three cation groups 
having different bonding geometries; (i) three Mg atoms, 
(ii) nine Zn atoms that are bonded to Oivig, and (iii) six 
Zn atoms that arc not bonded to OMg- For each cation 
group, the four nearest anion distances are slightly differ- 
ent, (i) ax'+3eq, (ii) ax+2eq+eq' , and (iii) ax+3eq, using 
the notations in Fig.fTUl However, the groups (i) and (iii) 
have same bond angles about the cation. It might be pos- 
sible that the Zn atoms (ii and iii) have non-uniform co- 
ordination geometries depending on the proximity to the 
Mg, but it is hardly expected that Zn (iii) and Mg (i) have 
the same bond angles. In this regard, the coordination 
geometries of Zn and Mg are not sensibly distinguished 
in the anion-lattice deformation model. In Fig. [121 the 
Zn— O and Mg— O bond distances in Zno.85Mgo.15O are 
plotted as obtained from the cation-lattice deformation 
model. The refined structural parameters (a, c, u) were 
used to calculate the bond distances. For both Zn and 
Mg, the first coordination shells consist of one longer 
(c- axial) bond and three shorter (equatorial) bonds. As 
reflected by the difference between the longer and shorter 
bonds, Mg atoms are found to sit in more regular tetra- 
hedra of oxygen than are Zn atoms. 



IV. DISCUSSION 

For various AX-type compounds including ZnO and 
MgO, the relative stabilities of wurtzite, zinc blende, and 
rock salt structures have been a long-standing and in- 
triguing subject The correlations between the bi- 
nary composition and the favored structure type have 
been proposed in several different ways, but commonly 
using the bond ionicity and the ion sizes of A and X as 
the principal parameters . ^^'^^1^^ If we limit our attention 
to ZnO and MgO, the observed crystal structures can be 
explained simply by the ionicity consideration, in other 
words, electronegativities of Zn and Mg. With interme- 
diate electronegativity, Zn can achieve the tetrahedrally 
directed valence orbital via facile sp'^-hybridization. On 
the other hand, Mg is more electropositive and its va- 
lence orbital is dominated equally by the 3px, 3py, and 
3pz components, by which the nearest oxide ions are ori- 
ented in the octahedral geometry. Hence Zn and Mg 
adopt the four- and six-fold coordinations, respectively, 
in their binary oxides. Given those coordination prefer- 
ences, the alloying with MgO should destabilize the host 
ZnO lattice. Indeed the thermodynamic solubility limit 
in Zni_a;Mg2;0 is reached at a low Mg-concentration, 

«15%.2'44 

The ZnO has wurtzite parameters of c/a = 1.6021 
and u = 0.3829, which are ideally 1.633 and 0.375, 
respectively.'^ Considering the range of c/a (1.600^1.645) 
observed from the binary wurtzites;^ ZnO is near the 
borderline of wurtzite stability. On the other hand, the 
wurtzite MgO structure is not available experimentally. 
Instead, the first-principle methods based on density 
functional theory (DFT) have predicted that the wurtzite 
MgO will have much smaller c/a than those of existing 
wurtzite phases. An earlier DFT study has reported an 
extremely distorted hexagonal structure {c/a = 1.20, u 
= 0.5) for the wurtzite MgO, in which Mg sits on a mir- 
ror plane and has a bi-pyramidal coordination.*^ How- 
ever recent studies, independently conducted by Janotti 
eiaL*^ and Gopal and Spaldin^^ have reported c/a = 
1.51 and u = 0.398. More relevant to our experimen- 
tal work, Malashevich and Vanderbilt performed the ge- 
ometry optimization of Zni_a;Mga;0 supercells {x = 0, 
i, J, i, and 5)'^ They showed that the c/a decreases 
with increase of Mg content, and u(Mg) is larger than 
u(Zn) in all the supercell compositions. Therefore those 
DFT studies on MgO and Zni_3;Mga;0 consistently sur- 
mise that the Mg-substitution in ZnO would result in the 
decrease of c/a, in agreement with the previously estab- 
lished idea that the c/a deviates farther from the ideal, 
when the bonding character becomes more ionici^ Our 
Rietveld study on Zni_j;Mga;0 has also shown that c/a 
is gradually decreased with the Mg-substitution.^ The 
Mg-substitution in ZnO is not disfavored from packing 
considerations, but may reduce the electrostatic stabi- 
lizations in the long range through the change in c/a 
ratio. At the solubility limit Zno.85Mgo.15O, the c/a is as 
low as 1.600, and this value seems to represent the lowest 
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c/a that can be sustained by the wurtzite lattice. 

The energy-minimized Zni_a;Mga,0 supercells of Mala- 
shevich and Vanderbilt^ have w(Mg) in the range 
0.387'-0.390, and the u(Zn), 0.378-0.382, which is not 
in agreement with the PDF results shown in Figs. [9] and 
fm As a consequence of u(Mg) > u(Zn), the sponta- 
neous polarization in the computed Zni_a;Mg2;0 super- 
cells become larger for the compositions with higher Mg 
content)^ while the Rietveld analysis indicated that the 
u parameter of Zni_2;Mg2;0 average structures and there- 
fore the spontaneous polarization in the Clausius-Mosotti 
limit decrease upon Mg-substitution.^ The above discrep- 
ancies are in fact connected with a more fundamental 
question of how the tetrahedral geometries of Zn04 and 
Mg04 would differ in the extended solids. Understand- 
ing of cation local geometries in Zni_2:Mg2:0 is impor- 
tant also for the application of ZnO/Zni_a;Mga;0 het- 
erojunction devices. While tetrahedral Mg04 is found 
in a few minerals such as MgAl204, Ca2MgSi207)^ and 
K2MgSi50i2r^ Zni_^Mg^O provides an opportunity to 
explore its equilibrated tetrahedral geometry. However 
the computational and experimental results on the inter- 
nal tetrahedral geometry of Zni_a;Mga;0 do not agree 
well. In the supercell structures of Malashevich and 
Vanderbilt,— the Zn04 tetrahedra are more regular than 
the Mg04 tetrahedra (as judged from the differences of 
the cation— anion bond distances), in contrast to the PDF 
findings presented in Fig. [121 These points may be fur- 
ther clarified by using alternative techniques such as neu- 
tron PDF or extended x-ray absorption fine structure. It 
is equally possible that the DFT methods must be re- 
examined, with inclusion of perhaps more accurate ex- 
change functionals. 

V. CONCLUSION 

Raman and synchrotron x-ray PDF have been em- 
ployed to probe the local structures of Zn and Mg atoms 



that occupy the common crystallographic site in the 
wurtzite alloys Zni-j^Mgj^O (x — 0, 0.05, 0.10, and 
0.15). Regardless of the Mg concentration, structure 
model, and r-range, the PDF refinements consistently 
show that the Mg atoms have smaller out-of-center tetra- 
hedral displacements than the Zn atoms. Even for the 
Zno.95Mgo.05O that has only 5 mol% MgO, the atomic co- 
ordinates of Zn and Mg were similarly obtained as those 
for Zno.9oMgo.ioO and Zno.85Mgo.15O, thereby demon- 
strating the fidelity of PDF technique. It is understood 
that the lessened tetrahedral distortion of Mg04, com- 
pared with Zn04, leads to the decrease of wurtzite u 
parameters in the Zni-j^Mg^jO average structures. The 
hexagonal c/a ratio of Zni_a;Mg3:0 decreases with the 
increase of Mg concentration, which is ascribed to the 
ionic character contributed from MgO. Therefore, from 
both ionicity and dipole strength viewpoints, we expect 
that the Zni-^iMg-cO alloys will have smaller sponta- 
neous polarizations than ZnO. It is noted that in terms of 
Zn04 and Mg04 local geometries, our experimental re- 
sults are somewhat contrary to the descriptions obtained 
from DFT-based studies. 
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